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EXPERIMENTAL AND THEORETICAL RESULTS FOR BUCKLING 
OF ECCENTRICALLY STIFFENED CYLINDERS 1 

By Michael F. Card and Robert M. Jones 
Langley Research Center 

SUMMARY 

Results from an experimental and theoretical study of the effect of stiffener eccen- 
tricity (one-sidedness) on buckling are reported. In the experimental investigation, 
axial -compression tests were conducted on twelve longitudinally stiffened cylinders which 
represent six configurations with internal or external, integral or Z-stiffeners. For cer- 
tain configurations externally stiffened cylinders were found to carry over twice the load 
sustained by their internally stiffened counterparts. A solution was obtained by a 
Galerkin procedure for buckling of clamped, eccentrically stiffened, orthotropic cylinders 
under axial compression, lateral pressure, or hydrostatic pressure. The experimental 
results for axially loaded cylinders range from 70 to 95 percent of the corresponding 
theoretical predictions. However, it is concluded from the comparison that a membrane, 
prebuckled shape is sufficient for assessment of eccentricity effects. Furthermore, com- 
parison of solutions for clamped and simply supported cylinders with the test data and 
with previous theoretical results reveals that edge clamping has a significant effect. 

INTRODUCTION 

The influence of eccentricity or one-sidedness of stiffeners on the buckling strength 
of stiffened cylindrical shells has been demonstrated in several studies (refs. 1 to 6). 
Buckling predictions derived from these studies indicate that the magnitude of the eccen- 
tricity effects can be very substantial for certain configurations. Unfortunately, few 
experiments have been devised to verify the trends of those theories. A need conse- 
quently exists for buckling test data on eccentrically stiffened cylinders. Furthermore, 
correlation of test data with available theory is limited because the theoretical boundary 
conditions usually considered are simple support whereas the test cylinders may have 
nearly clamped edges. - 

1 A part of the material presented herein is contained in the proceedings of 
AIAA/ASME Seventh Structures and Materials Conference (Cocoa Beach, Florida), 

April 18-20, 1966, pp. 23-34, under the title "Buckling of Axially Compressed Cylinders 
With Eccentric Longitudinal Stiffeners," by Michael F. Card and Robert M. Jones. 



The purpose of the present paper is 
to report the results of an experimental 
and theoretical study of the effects of stiff- 
ener eccentricity on the buckling strength 
of stiffened cylinders. (See fig. 1 for con- 
figuration.) In the experimental investiga- 
tion, axial-compression tests were con- 
ducted on twelve longitudinally stiffened 
cylinders which represent six configura- 
tions with internal or external integral or 
Z- stiffeners. Preliminary results of the 
experimental investigation are reported in 
reference 7. 

In the present paper, the test results 
are compared with theoretical predictions 
for simply supported and clamped cylinders in order to assess the effect of boundary con- 
ditions on the buckling load. The theoretical predictions are based on the buckling theory 
given in reference 6. The theory presented therein is a generalization of that given in 
reference 2. The solution for simply supported edges is reported in reference 6 and 
the solution for clamped edges obtained by use of a Galerkin procedure is reported in 
appendix A of this paper. The solution in appendix A is valid for clamped, eccentrically 
stiffened, orthotropic cylinders under axial compression, lateral pressure, or hydrostatic 
pressure. Some additional results on the influence of boundary conditions are presented 
in appendix B. 



Figure 1.- Stiffened-cylinder configuration. 


SYMBOLS 

The units for the physical quantities defined in this paper are given both in U.S. 
Customary Units and in the International System of Units (SI) (ref. 8). Factors relating 
these two systems are presented in appendix C. 

a ring spacing (see fig. 1) 

A cross-sectional area of stiffener 

At total cross-sectional area of stiffeners and cylinder 

b stringer spacing (see fig. 1) 

Ef3 

D bending stiffness of isotropic plate, — 7-= -r 

12 (l - p 2 j 
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D x ,Dy 


D xy 

E 


ExjEy 


G 

G xy 

h 

I 


J 


kx 

L 

m 

M x 

n 


N x 

Nx>Ny 


N e 

N c 


bending stiffnesses of orthotropic cylinder in longitudinal and circumferen- 
tial directions, respectively 

twisting stiffness of orthotropic cylinder 

Young's modulus 

extensional stiffnesses of orthotropic cylinder in longitudinal and circumferen- 
tial directions, respectively 

E 

shear modulus, — 7 

2(1 + p) 

in-plane shear stiffness of orthotropic cylinder 

height of stiffener element 

moment of inertia about centroid of stiffener 


torsional constant for stiffener, V tjL h(i - ®li§e j an j 1 _|2L) ( see re f. 9) 

Lj \3 ir & h 2t se l 


axial-compression stress coefficient, 


cylinder length 


a tl/ 
Dtt 2 


number of longitudinal buckle half-waves 


additional moment per unit length during buckling 


number of circumferential buckle waves 


experimental value of n 


additional axial force per unit length of circumference during buckling 
applied axial and circumferential forces per unit length, respectively 


experimental value of N x 


theoretical value of N x for clamped edges 
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Ns theoretical value of N x for simply supported edges 

p c corrugation width (see fig. 2) 

P e experimental value of axial buckling force 

R radius of cylinder middle surface 

t cylinder thickness 

A A 

t effective thickness of stiffened isotropic cylinder, £ + t or + t 

a b 

t c corrugation thickness (see fig. 2) 

t se thickness of stiffener element 

u,v,w displacements from membrane prebuckled shape in longitudinal, circumferen- 
tial, and radial directions, respectively 

x,y,z longitudinal, circumferential, and radial coordinates, respectively 

T 2 

Z curvature parameter, — 

Rt 

z distance from stiffener centroid to cylinder middle surface (see fig. 1), 

positive when stiffener on outside 

6 axial strain 

cr axial stress 



0g experimental average axial buckling stress 

6>c corrugation angle (see fig. 2) 

p. Poisson's ratio 

p x ,p y Poisson's ratios for bending of orthotropic cylinder in longitudinal and circum- 

ferential directions, respectively 
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Poisson's ratios for extension of orthotropic cylinder in longitudinal and 
circumferential directions, respectively 

Subscripts: 

r,s denote properties of rings (circumferential stiffening) and stringers 

(longitudinal stiffening), respectively 

When the subscripts x and y follow a comma, they indicate partial differentia- 
tion of the principal symbol with respect to x and y. 

EXPERIMENTAL INVESTIGATION OF LONGITUDINALLY 
STIFFENED CYLINDERS 

Specimens 

Buckling data and stress- strain curves were obtained for twelve longitudinally 
stiffened cylinders. Four specimens of about 9.55-inch (24.26-cm) radius and about 
0.0279-inch (0.0709-cm) thickness employed integral rectangular stiffeners; eight 
specimens of about 15.92-inch (40.44-cm) radius and about 0.0399-inch (0.1013-cm) 
thickness employed Z- stiffeners which were riveted to the cylinder. For each configura- 
tion, two identical cylinders were constructed, one with external stiffening and the other 
with internal stiffening. All specimens were designed with closely spaced stiffeners to 
preclude local buckling of the cylinder. Construction details are given in figure 2 and 
table I. 

The cylinders with integral rectangular stiffeners were machined from 3/8-inch 
(0.95-cm) aluminum plates. The dimensions given in table I were obtained by careful 
measurement. The cylinder thickness t is the average of a large number of microm- 
eter measurements. The scatter in individual thicknesses was ±8 percent of the tabulated 
value and is attributed to difficulties in measurement on the specimen surface. Note in 
table I, however, that the average-thickness values are fairly uniform. There were two 
butt joints per cylinder, each with two doubler plates 0.70 inch (1.78 cm) wide and 
0.020 inch (0.064 cm) thick fastened to the cylinder with a single row of 3/32-inch 
(0.24-cm) rivets spaced at 5/16 inch (0.79 cm) on each side of the joint. 

The Z- stiffened cylinders were nominally identical to certain test specimens dis- 
cussed in reference 10 except that the Z-stiffeners were located on the external surface of 
the present specimens, whereas in reference 10 the Z-stiffeners were located on the 
internal surface. Four externally stiffened cylinders were tested to complement the 
results of the four comparable internally stiffened cylinders reported in reference 10. 

All dimensions are nominal except the area of the stiffened cylinder cross section A^, 
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(a) Riveted Z-stiffened cylinder. 


R = 9l55 
(24.3) 

STIFFENER SPACING = 1.00(2.54) 


(b) Integral rectangular stiff- 
ened cylinder. 
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which was determined by weighing and the 
cylinder thickness t, which was obtained 
by micrometer measurements. Individual 
cylinder thicknesses were within the usual 
tolerances specified for rolled sheet. 
There were three lap joints 3/4 inch 
(1.91 cm) wide, each with a double row 
of 3/32-inch (0.24-cm) rivets spaced at 
5/ 16-inch (0.79-cm) intervals. The 
rivet spacing for the stringers was also 
5/16 inch (0.79 cm). 

Procedure 


(c) Ring- stiffened, corrugated cylinder. 

Figure 2.- Dimensions of stiffened cylinder cross sections. All dimen- 
sions are in inches (centimeters). 


The cylinders were loaded in axial 
compression in the 1,200-kip-capacity 
(5.3-MN) universal static testing machine 
at the Langley Research Center. The ends 
of the cylinders were ground flat and 
parallel after being potted in a soft, low-temperature alloy for ease in machining. The 
cylinders were carefully alined in the testing machine to insure uniform bearing between 
the ends of the cylinders and the platens of the testing machine. Circularity of the cylin- 
ders was maintained during the tests by the use of 1/4-inch-thick (0.64-cm) aluminum 
bulkheads which were held in place by the potting material at each end of the cylinder. 


Twelve resistance-wire strain gages with a 6-inch (15-cm) gage length were 
mounted in back-to-back pairs on the cylinder and on the stringers to determine stress 
distribution and to detect buckling. The overall shortening of the distance between the 
testing machine platens was measured with the use of strain gages mounted on a small 
cantilever beam, the fixed end of which was fastened to the lower platen. A long rod con- 
nected the beam to the upper platen shown in figure 3(a) so that the shortening of the 
cylinder was equal to the deflection of the cantilever. The strain and shortening measure- 
ments were recorded during the test at a virtually continuous rate on the Langley central 
digital data recording system. 


Stress Distribution 

Stress- strain curves derived from the strain-gage data for each of the test cylin- 
ders are presented for stresses below the buckling stress in figure 4. The strain plotted 
is that measured at the location which experienced the maximum strain during the test. 
The value of the strain corresponds to the average of data from a pair of back-to-back 
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TABLE I.- DIMENSIONS AND RESULTS 


(a) U.S. Customary Units 


Cylinder 

Aluminum 

Type of 

Stiffener 

t, 

L, 

R, 

At, 

p e , 


n e 

N e , 

N c , 

N s , 

alloy 

stiffener 

location 

in. 

in. 

in. 

in 2 - 

kips 

ksi 

kips/in. 

kips/in. 

kips/in 

1 

2024-T351* 

Integral 

External 

0.0283 

38.00 

9.55 

3.69 

112.6 

30.5 

6 

1877 

2036 

1145 

2 

2024-T351* 

Integral 

Internal 

.0277 

38.00 

9.55 

3.72 

48.0 

12.9 

6 

800 

997 

687 

3 

2024-T351* 

Integral 

External 

.0275 

23.75 

9.55 

3.70 

127.2 

34.4 

7 

2120 

2595 

1326 

4 

2024-T351* 

Integral 

Internal 

.0280 

23.75 

9.55 

3.63 

58.0 

16.0 

6 

965 

1226 

722 

5 

7075-T6* 

Z 

External 

.0410 

59.00 

15.92 

7.15 

169.6 

23.7 

5 

1696 

2289 

1193 

6** 

7075-T6* 

Z 

Internal 

.0401 

59.00 

15.92 

7.25 

99.5 

13.7 

6 

995 

1174 

766 

7 

7075-T6* 

z 

External 

.0391 

41.83 

15.92 

7.02 

186.3 

26.6 

7 

1863 

2617 

1283 

8** 

7075-T6* 

z 

Internal 

.0401 

41.83 

15.92 

7.31 

126.0 

17.2 

6 

1260 

1383 

773 

9 

7075- T6* 

z 

External 

.0392 

29.68 

15.92 

7.06 

273.0 

38.7 

7 

2730 

3426 

1574 

io*+ 

707 5- T6* 

z 

Internal 

.0403 

29.68 

15.92 

7.21 

171.0 

23.7 

7 

1710 

1813 

838 

11 

7075-T6* 

z 

External 

.0395 

20.77 

15.92 

7.18 

378.0 

52.6 

8 

3780 

4947 

2068 

12** 

707 5- T6* 

z 

Internal 

.0401 

20.77 

15.92 

7.18 

273.0 

38.1 

8 

2730 

2845 

1027 


*E = 10.5 X 10 3 ksi; n = 0.32. 
**Data are taken from reference 10. 


(b) International System of Units 


Cylinder 

Aluminum 

Type of 

Stiffener 

t. 

L, 

R, 

At? 

p e , 

^e> 


N e , 

N c , 

N s , 

alloy 

stiffener 

location 

cm 

cm 

cm 

cm 2 

kN 

MN/m 2 

u e 

MN/m 

MN/m 

MN/m 

1 

2024-T351* 

Integral 

External 

0.0719 

96.5 

24.26 

23.8 

500 

210 

6 

329 

357 

200 

2 

2024-T351* 

Integral 

Internal 

.0704 

96.5 

24.26 

24.0 

214 

89 

6 

140 

175 

120 

3 

2024-T351* 

Integral 

External 

.0699 

60.3 

24.26 

23.9 

566 

237 

7 

371 

454 

232 

4 

2024-T351* 

Integral 

Internal 

.0711 

60.3 

24.26 

23.4 

258 

110 

6 

169 

215 

127 

5 

707 5- T6* 

Z 

External 

.1041 

149.9 

40.44 

46.1 

754 

163 

5 

297 

401 

209 

6** 

707 5- T6* 

Z 

Internal 

.1019 

149.9 

40.44 

46.8 

443 

94 

6 

174 

206 

134 

7 

707 5- T6* 

Z 

External 

.0993 

106.2 

40.44 

45.3 

829 

183 

7 

326 

458 

225 

8 ** 

707 5- T6* 

z 

Internal 

.1019 

106.2 

40.44 

47.2 

560 

119 

6 

221 

242 

135 

9 

707 5- T6* 

z 

External 

.0996 

75.4 

40.44 

45.5 

1214 

267 

7 

478 

600 

276 

io** 

707 5- T6* 

z 

Internal 

.1024 

75.4 

40.44 

46.5 

761 

163 

7 

299 

317 

147 

11 

707 5- T6* 

z 

External 

.1003 

52.8 

40.44 

46.3 

1681 

363 

8 

662 

866 

362 

12** 

707 5- T6* 

z 

Internal 

.1019 

52.8 j 

40.44 

46.3 

1214 

263 

8 

* 478 

498 

180 


*E = 72 400 MN/m 2 ; fx = 0.32. 
**Data are taken from reference 10. 
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(a) Integrally stiffened cylinder no. 3. 



gages mounted on the stiffeners or cylinder. 

The strain distribution was nearly uniform for 
all tests. Individual strain-gage readings at 
loads near buckling differed at most by 10 per- 
cent of the plotted values. Strains measured in 
back-to-back gages prior to buckling differed by 

about 1^ percent of their average value. 

& 

All the stress-strain curves of figure 4 
are linear up to at least 85 percent of the maxi- 
mum stress carried by the cylinder. The slopes 
of the linear portion of the curves are in reason- 
able agreement with conventional values of 
Young's modulus. Slopes of the load- shortening 
curves derived from the average of the canti- 
lever measurements were also in agreement 
with accepted modulus values. The strain data 
indicate no local buckling of the cylinders prior 
to general buckling. The data also indicate that, 
for all cylinders except cylinder 4, the buckling 
load and maximum load could be considered 
identical. For cylinder 4, strain reversal 
occurred at 94 percent of the maximum load and 
was considered to determine the buckling load. 

Buckling 


(b) Integrally stiffened cylinder no. 1.. 



(c) Z-stiffened cylinder no. 11. 

Figure 3.- Buckled cylinders. L-66-4432 


For the externally stiffened cylinders, 
buckling was accompanied by a loud report and 
the appearance of several diamond- shaped 
buckles of large amplitude. In contrast, for the 
internally stiffened cylinders, buckling was a 
milder phenomenon with similar buckles of 
smaller amplitude. Photographs of buckled 
cylinders are shown in figure 3. The diamond- 
shaped buckles occurred in only one tier in all 
cases except cylinder 1 for which two tiers 
developed. (Compare figs. 3(a) and 3(b).) The 
number of circumferential waves % into which 
the cylinders buckled experimentally is shown 
in table I. 
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(a) Integrally stiffened cylinders. 



(b) Z-stiffened cylinders. 

Figure 4.- Stress-strain curves. 

The buckling loads together with the corresponding average axial compressive 
stresses carried by the cylinder at buckling are given in table I. Data for cylinders 6, 8, 
10, and 12 from reference 10 are included for comparative purposes. Examination of the 
buckling loads reveals that, for the test boundary conditions, locating the longitudinal 
stiffening on the external surface of the cylinder can produce a significant increase (to 
over twice as much) in axial load-carrying ability over that of the comparable internally 
stiffened cylinder. 

THEORETICAL RESULTS AND COMPARISON OF THEORY AND EXPERIMENT 
FOR LONGITUDINALLY STIFFENED CYLINDERS 

The influence of the boundary conditions of the present longitudinally stiffened 
cylinders must be considered because the cylinders are effectively short as will be seen 
subsequently herein. For this reason, the experimental results are compared with 
buckling predictions from solutions for cylinders which have, in one case, simply sup- 
ported edges and, in the other, clamped edges. The calculations for simply supported 
cylinders were obtained by employing equation (15) of reference 6, the boundary condi- 
tions of which can be identified as N x = v = w = M x = 0. The calculations for clamped 


9 




cylinders (u = v = w= w x = 0) were obtained from the solution in appendix A. Note 
especially the difference in axial constraint between the two sets of boundary conditions. 

The solution for clamped cylinders presented in appendix A is based on the stabil- 
ity differential equations of reference 6 wherein a membrane prebuckled state was uti- 
lized. The buckling displacements, u, v, and w, were expanded in trigonometric 
series, each term of which satisfies the clamped- edge boundary conditions. A general- 
ized Galerkin procedure was applied to obtain a set of algebraic equations which only 
admit solutions for certain values (eigenvalues) of the load parameter. Buckling loads 
and mode shapes correspond to the lowest eigenvalue and its eigenvector, respectively. 


Theoretical and experimental results for integrally stiffened and Z- stiffened 
cylinders are shown in figures 5 and 6, respectively, and in table I. The theoretical 
results in figures 5 and 6 are based on the average of the stiffened- cylinder dimensions 
given in table I. In figures 5 and 6, substantial boundary- condition effects are evident. 
Note that the percentage difference between the solutions for simply supported and 

clamped edges is greater for 
externally stiffened cylinders than 
for their internally stiffened 
counterparts. Thus, boundary 
conditions are more significant 
for externally than for internally 
stiffened cylinders of the configura- 
tions considered. This result 
should not be interpreted too gen- 
erally because the opposite effect 
is shown in appendix B for certain 

BUCk t^in^rraS al r tiffened Cylinders ring-stiffened corrugated cylinders. 



Figure 5.- 
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Figure 6.- Buckling of axially compressed Z- stiffened cylinders (as 
described in table I and fig. 2). 


The theoretical buckling 
results shown in figures 5 and 6 
have been extended to large values 
of the length- radius ratio L/R to 
investigate the convergence of the 
solution for clamped edges to the 
solution for simply supported 
edges. Note from the figures that 
even for L/R = 10 substantial 
boundary- condition effects are 
present for externally stiffened 
cylinders. Although it is expected 
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that boundary conditions have an influence in stiffened cylinders, the actual magnitude of 
the effects at such large values of L/R is somewhat unexpected. In particular, 
boundary- condition effects are present in longitudinally stiffened cylinders at values of 
L/R for which the boundary conditions of unstiffened cylinders have usually been 
regarded as immaterial. Appendix B shows that the presence of these effects is due in 
part to the additional axial constraint implied by the present clamped boundary conditions. 
However, appendix B also indicates that other stiffened configurations have large 
boundary- condition effects at relatively large values of L/R. 

In the course of the theoretical investigation, it was observed that the torsional 
stiffness G S J S had a much larger influence on the buckling load for integrally stiffened 
cylinders than for Z- stiffened cylinders. This situation is explained by noting that the 
integrally stiffened cylinders have a relatively shallow stiffener so that the ratio of 
twisting to bending stiffness G S J S /E S I S is much larger than for the Z- stiffened cylin- 
ders. To assess the magnitude of this effect, calculations were made in which G S J S 
was set equal to zero. Whereas the resulting buckling loads for Z- stiffened cylinders 
were only 1 percent lower than the loads shown in figure 6, the corresponding values for 
integrally stiffened cylinders were 10 to 20 percent lower than those shown in figure 5. 
Thus, for some types of stiffening, it is possible to arrive at erroneous conclusions 
regarding the behavior of stiffened shells if torsional stiffness is neglected. 

In general, the qualitative trends of the theory are verified by the experimental 
results. Evidently, the test cylinders were more nearly clamped than simply supported. 
Quantitatively, the experimental data vary from 70 to 95 percent of the theoretical pre- 
dictions for clamped cylinders. Better agreement occurs for internally than for exter- 
nally Z- stiffened cylinders and for externally than for internally integrally stiffened 
cylinders. In addition to boundary conditions, contributions to the difference between 
theory and experiment can result from neglect of initial imperfections and from the 
absence of an accurate prebuckling solution in the theory. However, comparison of 
theoretical and experimental results suggests that use of a membrane prebuckling shape 
is sufficient for assessment of eccentricity effects. 

CONCLUDING REMARKS 

The results of an experimental and theoretical study of the effects of stiffener 
eccentricity on the buckling strength of stiffened cylinders have been discussed. The 
buckling experiments were conducted on 12 axially compressed, longitudinally stiffened 
cylinders in order to study the relative effect of locating stiffeners on the internal or 
external surface of the cylinder. Externally stiffened cylinders were found to carry axial 
loads up to twice those sustained by their internally stiffened counterparts, a fact which 
has previously been predicted on theoretical grounds. 
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The experimental buckling data for axial compression of longitudinally stiffened 
cylinders have been compared with an existing theoretical solution for simply supported 
cylinders and with a solution described in an appendix for clamped stiffened cylinders. 
Quantitatively, the experimental data vary from 70 to 95 percent of the theoretical pre- 
dictions for clamped cylinders. The comparison indicates that use of a membrane pre- 
buckled shape is sufficient for assessment of eccentricity effects in longitudinally stiff- 
ened cylinders. It was also found that for some configurations, the torsional stiffness of 
stiffeners can have a substantial effect on the buckling load. In addition, large effects of 
edge clamping and axial constraint were found for stiffened cylinders with length- radius 
ratios for which unstiffened cylinders are generally considered to be insensitive to 
boundary conditions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 7, 1966, 

124-11-06-04-23. 


12 



APPENDIX A 


BUCKLING OF CLAMPED ECCENTRICALLY STIFFENED ORTHOTROPIC 
CYLINDERS UNDER AXIAL COMPRESSION, LATERAL PRESSURE, 
OR HYDROSTATIC PRESSURE 


The following stability equations for eccentrically stiffened orthotropic cylinders 
which have the configuration shown in figure 1 are presented in reference 6. The equa- 
tions are expressed in terms of displacements u, v, and w from a membrane pre- 
buckled shape which has externally applied force resultants N x and Ny: 


Ex 


E s A s 


t 1 yEx 


^y E J 


+ “b-j u ’ xx + ***>” + [** + + R(1 . ^ 


_ EqAq 

“ z s — H w ,xxx = 0 
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Ey E r A r \ ( ^x^v \ 

.1 - tw + —h? * G *y V 'XX + K + I^d”,xy 
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E v E r A r \ 
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_ E r A r 

Z r — W } yyy = 0 
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r^r 
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Mx^y a 

EpAp 


_ E r A r 
"a >j 


+ — — I v y - z r v^yyy + 


Ey E r A r \ 

— + — ^ — w 


R2 U - a 


( — 2 — E r A r \ — 

\% - R z r — a~ J w ,yy + N xW,xx = 0 


(A3) 


The clamped- edge boundary conditions at x = 0 and L chosen for the present investiga- 
tion are 

u = v = w = w jX = 0 (A4) 
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APPENDIX A 


A solution which is apparently equivalent to solution of equations (Al) to (A3) with 
boundary conditions expressed by equation (A4) is obtained by an indirect method in ref- 
erence 3. For the sake of clarity and completeness, a straightforward solution is 
obtained here by a generalized Galerkin procedure. (See ref. 11.) The displacements 
are expanded in a series of functions which satisfy the boundary conditions, equation (A4): 



The Galerkin procedure requires that the following equations be satisfied: 


n 2 ttR , 

- L * ^ 

) 1 
Jq 

l UiQx(u,v,w)dx dy = 0 

'o 

p 2irR , 

- L 

J o 

) Vi Qy(u,v,w)dx dy = 0\ 

n 2 vR , 

” L . 

) 1 
; 0 

\ WiQ z (u,v,w)dx dy = 0 

’o 

J 


(A6) 


where Qx, Qy, and Q z are the left-hand sides of equations (Al), (A2), and (A3) with 
the series, equations (A5), substituted for u, v, and w, respectively. If I terms are 
taken in each of the series in equations (A5), then equations (A6) represent a set of 31 
linear, homogeneous, algebraic equations in 31 unknowns. The set is symmetric if the 
reciprocal relations 

fJ-xEy = / i y E x 


fi x D y = fi y D x 


are employed. Examination of the resulting equations reveals that the equations with 
a^ (i odd), bj (i even), and Cj (i even) are not coupled with those which involve 
aj (i even), bj (i odd), and c^ (i odd). Thus, there are two independent sets of 
simultaneous equations, one of which represents a symmetric w-deformation (i odd 
in cjJ and the other an antisymmetric w-deformation (i even in cj). Both sets can be 
written as 

AX = BX ( A ' 
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1 

a 22 

a 33 

• 

• 

B 11 b 12 b 13 • • • 
b 21 b 22 b 23 • • • 

b 31 b 32 b 33 • • • 

• • • 

• ♦ * 

• • • 

Cll c 12 

C21 c 22 C23 

C32 C33 C34 

• • • 

• « 

• # 


B 11 

d 22 

d 33 

• 

• 

En e 12 e 13 • • 

e 21 e 22 e 23 • • 

e 31 e 32 e 33 • • 

* 

• • • 

• • • 

(symmetric about 
main diagonal) 


F 11 f 12 
f 21 f 22 f 23 

f 32 f 33 f 34 


(zero) 


(zero) 


(zero) 


(symmetric about 
main diagonal) 


(zero) 


(zero) 


P 11 p 12 

p 22 p 23 

p 33 P 34 
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and 


X 


symmetric = 


a 2 


a l 

a 4 


a 3 

a 6 

• 


a 5 

• 

• 

• 

bl 


• 

b 2 

b 3 


b 4 

b 5 

• 

•^antisymmetric = 

b 6 

• 

ci 


c 2 

c 3 


c 4 

c 5 

• 


c 6 

• 

• 

• 


• 


(A10) 


For the symmetric case, the elements of the submatrices in equations (A8) and (A9) are 




tj _ 8 n i(2j - 1) 


/IjjEy 


li-'LR(2 1 )2.(2 i .i)2^1-MiMy + Gxy / 

e s a s /0 .s3/ it\ 3 




_ EeA 

Cij = z s 


Cij = 0 


~\ 

(i = j) 


(i = j)> (All) 
(i = j - 1) 


0 = 3 + 1), 


(Equations continued on next page) 
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f 


«< 


where 
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E\ 


E r A r \/n'2 


2/7T\ 2 


MxMy 


I ) + G * y (2i - « (0 


(i - i) 


Ei3=i(2i- 1) 


1 1 

7 % ] E r A r) 

n 

(2i - l) 2 - (2j - 2)2 (2i - l) 2 - (2j) 2 

\1 " MxMy a yi 

R 2 


+ z r 


E r A r /n \ 3 


a \R 


— I6A1 + 4A2 + 3Ag 

Fij = [(2j - 2) 4 + (2j) 4 ]Ai + [(2j - 2) 2 + (2j)^A 2 + 2A a 
Fij = -[(2j - 2) 4 Ai + (2j - 2) 2 A 2 + Agj 


(i = 3 * D 

(i = i - 1] 



(All) 
> Cone. 
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For the antisymmetric case, the elements of the submatrices in equations (A8) and (A9) 


are 


Am = 


E X Eg A, 


1] l 1 - MxMy 

8 n (2i - l)j 


/ MyE x ' 


(i = j) 




Cii- 0 


- E S A i= 

C-ij ~ -z s 




Bm = 


Ey ^ E r A r ^ n ^2 


1] l 1 - ' a 


F .. _ 8i 

El J " ir 


)(l) 2 + ^^^(l) 2 


(i = j) 
(i = j - 1) 
(i = j + 1) 

(i = j) 


(2i)2 - (2j - 1)2 (2i)2 - (2 j + l) 2 


> (A12) 


1 - 


E y + E r A r\ n 
M x My a JR 2 


+ - E r A r y3 


Jr a iR 


Fjj = [(2j - l) 4 + (2j + 1) 4 Jaj + [(2j - l) 2 + (2j + 1) 2 ]a 2 + 2A 3 (i = j) 
Fij = -[(2j - 1) 4 A X + (2j - 1) 2 A 2 + A 3 ] 

Pij = [(2J - l) 2 + (2J + D 2 ](f;) 2 Nx + 2 (|) 2 N y 


P ii = ' 


(2] - D 2 (l) 2nx + (r/^ 

where Aj, A 2 , and A 3 are the same as defined in equations (All). 


(i = i - 1) 
(i = j) 
(i = j - 1) 
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It is convenient for what follows to further abbreviate equation (A7) as 


A 

B 

c 


a 


0 

0 

(f 


a 


D 

E 

• 

b 

= 

0 

0 

0 

- 

b 

C T 

eT 

F 


c 


0 

0 

p 


c 


(A13) 


where the superscript T indicates the matrix transpose. Because of the zeros in the 
right-hand side of equation (A13), the following stability equation can be obtained by 
matrix operations: 


[f] - [c t |e t ] • 


"A 

|£1 

-1 

c 

_b t 

Idl. 


E_ 




(A14) 


where the superscript -1 indicates the matrix inverse. The left-hand side of equa- 
tion (A14) can be shown to be symmetric upon consideration of the symmetry properties 
of the coefficients A, B, C, etc. For specified combinations of the force resultants 
N x and Ny, [p] = x[p] where X is a scalar. For example, if the loading is axial 
compression, Ny = 0 and X = N x ; for lateral pressure, N x = 0 and X = p; and for 

hydrostatic pressure, N x = i Ny and X = p. In addition, if N x is fixed, the critical 
value of Ny can be found and conversely. In this manner, an interaction curve can be 
obtained if desired. 

The buckling load can be obtained from the stability equation (A14) by the following 
procedure. When the number of terms in the series expansions of the displacements as 
well as the number of buckle waves n in the circumferential direction is specified, 
the lowest eigenvalue can be extracted from equation (A14). The number n is varied 
over some range until an absolute minimum eigenvalue is found. In the preceding step, 
it is possible to obtain several relative minimums for the lowest eigenvalue; thus, solu- 
tions for a wide range of n values must be investigated before it can be assumed that 
the absolute minimum has been found. For a given problem, the buckling load is the 
lowest eigenvalue of the two absolute minimums obtained for the symmetric and anti- 
symmetric modes. 

Numerical values were obtained by use of a computer program. Convergence of 
the solutions was assured to an arbitrary accuracy by comparison of solutions with dif- 
ferent numbers of terms in the series expansions of the displacements. The number of 
terms necessary for convergence to an accuracy of three significant figures varied from 
three to eight as the shell length was increased. 
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INFLUENCE OF BOUNDARY CONDITIONS 


Buckling of Isotropic Cylinders 


The significance of the clamped- edge boundary conditions employed herein as 
applied to the buckling load (eq. (A4)) can be observed by comparing the present solution 
with exist ing solutions for axially loaded clamped cylinders and simply supported cylin- 
ders. The present solution can be specialized for unstiffened isotropic cylinders by 
taking the stiffener properties to be zero and making the following substitutions for the 
orthotropic material constants: 


Mx = My = M X = My = M 
Dx = Dy = Et 3 /l2 
Dxy = Gt 3 / 6 > 

E x = Ey C Et 
Qxy = Gt 


(Bl) 


The buckling solution thereby obtained is shown in figure 7 along with existing 
buckling solutions for a membrane prebuckled shape under various boundary conditions. 
Note that the solutions obtained in reference 12 for clamped and simply supported edges 



Figure 7.- Buckling of axially compressed isotropic 
cylinders. 


converge when Z is about 10. In contrast, 
the present results for clamped edges do not 
converge to the solution for simply supported 
edges until Z is about 100 because of the 
additional axial constraint implied by the 
present boundary conditions. The solution 
of reference 12 for simply supported cylin- 
ders can be shown to have the same boundary 
conditions as are used in the solution 
employed herein for simply supported stiff- 
ened cylinders (ref. 6). Thus, by analogy, 
it can be inferred that the substantial effect 
of clamped edge conditions on the buckling 
load of even relatively long stiffened cylin- 
ders shown in the present results is due in 
part to axial constraint. 
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Buckling of Axially Compressed Ring- Stiffened Corrugated Cylinders 

The present buckling equations can be specialized for ring- stiffened corrugated 
cylinders by substitution of the following stiffnesses in equations (Al) to (A3): 

Ex = 2Ejt c /(l + cos 6> c j] ^ 

Gxy = Gt c (Et c /Ex) 

Dx = Et c p c 2 jsin 2 6» c /3(l + cos 6> C )J 

% = D y = D xy = = My = M X = M y = 0, 

where t c is the corrugation thickness, p c is the corrugation width, and 9 C is the 
corrugation angle (t c , p c , and 6 C shown in fig. 2(c)). 

The effect on the buckling load of increasing cylinder length can be studied by 

obtaining the buckling load of cylinders the length of which is some multiple of the ring 

spacing. The results of such a procedure are 

shown in figure 8. In this case, the present 

clamped-edge solution converges to the simply 

supported- edge solution of reference 6 much more 

rapidly (on a percentage basis) for external than 

for internal rings. At a length- radius ratio L/R 

of 2, the clamped-edge solution is only 13 percent 

higher than the simply supported- edge solution for 

external rings; however, for internal rings, the 

corresponding value is 45 percent. Thus, boundary 

conditions are more significant for internal rings 

than external rings for corrugated cylinders of 

the configuration considered. This conclusion is 

interesting because, in contrast, the longitudinally 

stiffened cylinders considered were found to have 

substantially greater boundary effects for external Figure 8.- Buckling of axially compressed, ring- 

stiffened, corrugated cylinders (as described 

stiffeners than for internal stiffeners. in fig. 2). 


(B2) 
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APPENDIX C 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Con- 
ference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 8). 
Conversion factors for the units used herein are given in the following table: 


Physical 

quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 

Length 

in. 

0.0254 

meters (m) 

Force 

kips 

4.448 x 10 3 

newtons (N) 

Stress 

ksi = kips/in^ 

6.895 x 10 6 

newtons per square meter 


^Multiply value given in U.S. Customary Unit by conversion factor to 


{Wm 2 ) j 

obtain 


equivalent value in SI Unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

centi (c) 

10- 2 

kilo (k) 

10 3 

mega (M) 

10 6 
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